Universal persistence exponents in an extremally driven system 
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The local persistence R(t), defined as the proportion of the system still in its initial state at time 
t, is measured for the Bak-Sneppen model. For 1 and 2 dimensions, it is found that the decay of 
R(t) depends on one of two classes of initial configuration. For a subcritical initial state, R(t) ~ t~ e , 
where the persistence exponent 9 can be expressed in terms of a known universal exponent. Hence 
9 is universal. Conversely, starting from a supercritical state, R(t) decays by the anomalous form 
1 — R(t) ~ t TALL until a finite time to, where tall is also a known exponent. Finally, for the high 
dimensional model R(t) decays exponentially with a non-universal decay constant. 

PACS numbers: 05.70.Ln, 64.60.Ht 
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Extremally driven systems represent an important 
class of non-equilibrium models, with applications in- 
cluding sheared granular media, biological macroevolu- 
tion and interface depinning Their defining char- 

acteristic is that they are driven in the vicinity of an 
extremal quantity, often the minimum or maximum of an 
inhomogeneous scalar field. Such deterministic dynamics 
can be naturally realised in the limit of vanishing (ther- 
mal) noise Moreover, many extremal models are 
also 'critical' in an analogous sense to a continuous phase 
transition in equilibrium systems 

However, the temporal evolution of these systems is in 
many cases not fully understood. This is clearly demon- 
strated by recent work on one of the simplest, and there- 
fore most studied, extremally driven systems, namely the 
Bak-Sneppen model This has been shown to age 
in a similar manner to glassy systems |^,|l0|, suggesting 
that it does not reach a steady state in the timeframe of 
current simulation methods. This is in contrast to the 
usual claim that stationarity is reached 'after an exten- 
sive transient ' || . It now appears that this contradiction 
arises because most studies only measure functions of a 
single time variable, such as critical exponents etc., that 
can appear to become constant after a finite time even 
for a non-stationary system; a true test of stationarity 
(such as the aging result mentioned above) re quir e the 
measurement of two-time correlation functions |l l| , |l2| . 

The aim of this paper is to investigate one aspect of 
the temporal behaviour of the Bak-Sneppen that has not 
yet been studied, in the hope that this will further eluci- 
date the time evolution of extremally driven systems in 
general. The quantity that we consider is the local persis- 
tence R(t), defined as the proportion of the system at a 
time t that has not yet changed from its initial state. This 
quantity has been the focus of considerable attention in 
recent years, at least in part because it is often found 
to decay algebraically with a non-trivial persistence ex- 



ponent 9, R(t) ~ (const) t~ e , that (for equilibrium sys- 
tems) cannot be related to either the static or dynamic 
exponents [p~3| — |l6|l . The universality of 9 remains an open 
subject. In many cases it is not universal; that is, it may 
depend on the lattice coordination number or the precise 
choice of interaction term |l^,|l^]. However, it has been 
found to be universal in directed percolation fl6| . We 
find this to also be the case for the Bak-Sneppen model, 
and confirm our numerical finding by deriving an expres- 
sion for 9 in terms of the known universal exponent 7 [0 . 
We further find that, when starting from a supercriti- 
cal initial configuration, R{t) decays in an anomalous, 
non-algebraic manner. Finally, R(t) always decays ex- 
ponentially the infinite dimensional case. 

The model: The Bak-Sneppen model is defined as fol- 
lows A scalar quantity /, 6 [0, 1] is independently 
assigned to every site i of a (i-dimensional hypercubic 
lattice, which contains a total of N sites. For every 
unit time t, the smallest /j in the system is found, and 
it and its nearest neighbours are assigned new values, 
again drawn uniformly from [0,1]. Other initial condi- 
tions and interaction terms will be considered below. To 
generate the random values of the ft , we used L'Ecuyer's 
64 bit combined multiple recursive pseudo-random num- 
ber generator MRG63k3a, which has a quoted period 
of 2 377 The typical number of operations required 
to find the global minimum was reduced from O(N) to 
O(lnTV) by using a binary search tree, which also con- 
strains the system size to be a power of 2 |l9| . 

Typical results for a periodic ID lattice are presented 
in Fig. [j]. It can clearly be seen that the curves of R(t) 
for different system sizes collapse after rescaling the time 
scale by N, for all N > 2 10 sa 10 3 . For the largest sys- 
tem N = 2 21 , the simulation time was long and only one 
run was performed; however, for smaller systems R(t) 
was averaged over 10-100 different initial configurations. 
Also plotted in the figure is the theoretical prediction 
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R(t) oct~ 159 , which will be derived below. The numeri- 
cal R(t) appear to be converging towards this prediction 
at late times, but there is still a slight curvature to the 
lines (on log-log axes) even for the largest t we were able 
to simulate, making direct measurement of the exponent 
difficult. 

A similar picture emerges in two dimensions, as shown 
in Fig. ||. The theoretical prediction in this case R(t) oc 
t~ 2A3 (again, see below for details), which is consistent 
with the simulation data, although as before there is still 
some slight curvature to R(t) for the largest t available 
to our current computing resources. 

The universality of the exponent was tested by altering 
the microscopic details of the interactions. Three alter- 
native rules were considered : updating the second near- 
est neighbours as well as the nearest neighbours ("2nd 
NN's"); changing the value of the nearest neighbours 
from fi to \{fi + t ) , where the annealed random vari- 
able r is uniformly distributed on [0,1] ("Half NN"); and 
setting the minimum value to 1 but changing the near- 
est neighbours as normal ("Min to 1"). In all cases 9 
was found to be approximately the same. This has been 
tested in both one and two dimensions; the ID case is 
shown in Fig. 0. 

Theory: We shall now explain how the persistence expo- 
nent can be theoretically derived. This argument utilises 
known results for the Bak-Sneppen model, which for 
reasons of space will not be repeated here; the reader 
is instead referred to p|,|3|,pc| for their justification. It 
is known that the location of the minimum or 'active' 
site initially jumps around the lattice in an uncorrelated 
manner, but as the system evolves the lattice variables 
ft become spatially correlated, and the active site tends 
to remain in a localised region of the lattice for a finite 
amount of time before jumping to another, uncorrelated 
region of the lattice. Each such period of localised activ- 
ity is known as an 'avalanche.' 

Let the mean duration of an avalanche at time t be de- 
noted by (S(t)), and the mean number of sites changed 
at least once during the avalanche by (n cov (t)}, where 
the angled brackets '(■■■)' refer to averages over different 
initial configurations. By taking a suitably large system 
size N, these quantities will increase by an arbitrarily 
small amount between avalanches, allowing the contin- 
uum limit to be taken. It is known that both quantities 
increase algebraically in time 0, 



(S(t))~A 



(ncov(*)> ~ B 



t \ 7/(7-1) 



N 



N 



1/(7-1) 



(1) 

(2) 



Here, 7 is a universal exponent that depends on the lat- 
tice dimensionality and the symmetries of the interac- 
tions. We assume 7 > 1; The case of 7 = 1, realised in 
high dimensions, will be discussed later. 

By defintion, R(t) = 1 at t = and decreases by 
an amount 6R — 1/N whenever a lattice site changes 



its value for the first time. The average rate at which 
sites are changed by the avalanches is (n cov (t)) / (S(t)); 
however, only a fraction R(t) of sites have not already 
changed from their initial state. Recalling that the ac- 
tive site jumps to a random part of the system inbetween 
avalanches, then R(t) obeys the following equation 



dR(t) 



(n cov (t))R(t) 



dt (S(t))N 
Subsituting (£]) and (g) into @ gives 
R(t) ~ (const) r e 



(3) 



(4) 



with the persistence exponent 9 = B/A. 

Since A and B are non-universal constant prefactors 
to the algebraic growth of (S(t)) and (n cov (t)) respec- 
tively, one would expect that their ratio, and therefore 
9, is also non-universal; however, a peculiar feature of 
this model is that B/A is universal. This can be seen 
by combining together two exact equations, referred to 
as the gap and gamma equations in |^|, which gives the 
following equation relating (S(t)) and (n cav (t)), 



d(S(t)) _ (n cov (t)) 



dt 



N 



(5) 



Substituting (Q) and (§) into (§) gives the persistence 
exponent to be 



B 



A 7-1 



(6) 



Since 7 is universal, then so is 9. In 1 and 2 dimensions, 
7 = 2.70 and 1.70 respectively giving predicted val- 
ues of 8 ss 1.59 and 2.43, which were the values used in 
Figs. ||— - 0. As a further check on this analysis, note that 
d) and (§) predict that the quantity t{n cov (t))/(N{S(t)}) 
will approach the constant value B/A = 9 at late times. 
This is confirmed by the numerical results presented in 
Fig. g. 

The preceding argument breaks down when 7=1, 
which arises for dimensions above the upper critical di- 
mension al c , for which the most recent estimate is d c — 
4 plj . In this case, (S(t)) and (n cov (t)) both increase 
exponentially in time rather than algebraically [Q]. By 
repeating the same steps as above, it is straightforward 
to show that R(t) should now decay exponentially with 
a non-universal decay constant C, 



R(t) ~ (const) e" 



-ci 



(7) 



This is confirmed by the numerical results given in Fig. ^ 
for a system in which the minimum and another K — 1 
randomly chosen elements are assigned new values, which 
is the usual way of simulating an infinite dimensional sys- 
tem for this model. It is clear that the decay of R(t) is 
exponential with a decay constant that depends strongly 
on K. 
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Finally, we consider the effects of initial conditions. 
All of the preceding analysis and numerical results has 
assumed an initial configuration in which every /$ is 
uniformly distributed over the range [0,1]. This can 
be generalised by choosing an initial state in which the 
fi are uniformly distributed over the range [/o, 1] and 
are spatially uncorrelated. It is known that the Bak- 
Sneppen model has a 'critical value' f c £ (0, 1) such that 
the system is subcritical if /o < f c , and supercritical if 
fo > fc POl - (The value of f c is non universal and de- 
pends on the lattice structure and the microscopic inter- 
action term) . We have checked that any subcritical initial 
state gives rise to the same persistence exponent in both 
one and two dimensions, as predicted by our earlier argu- 
ment. However, for a supercritical initial configuration, 
the system enters into a single, infinite avalanche, and all 
of the preceding analysis fails. In fact, R(t) in this situa- 
tion is simply related to the growth of a single avalanche 
in time, which is already known to increase algebraically 
in time with a universal exponent tall 0- Thus R(t) 
does not decay algebraically but rather according to the 
expression 



1 - R(t) ~ (const) f 



(8) 



until a late time to when the avalanche 'touches' itself 
through the periodic boundary conditions. Confirmation 
of this scaling is given in Fig. | for one and two dimen- 
sions. In high dimensions, R(t) decays exponentially for 
both classes of initial state. 

Summary: We have found that the persistence R(t) in 
the Bak-Sneppen model can decay in one of three differ- 
ent ways. In low dimensions, it decays with a universal 
exponent 9 when starting from a subcritical state, but 
by the anomalous form 1 — R(t) cx t TALL when starting 
from a supercritical state. For high dimensions the decay 
is always exponential. We also provided a simple theo- 
retical prediction for the value of = 7/(7 — 1), which, 
when combined with the numerical results, gives strong 
confirmation of the universality of 9. 

It is not yet clear how this relates to the question of 
the existence of a steady state of the model, although 
it is interesting to note that 9 > 1 always, which is sig- 
nificantly higher than typical persistence exponents, and 
means that the number of untouched sites will rapidly de- 
cay to zero. Thus the observed non-stationarity must be 
due to temporal correlations within an avalanche, rather 
than the spread of avalanches throughout the system. 
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FIG. 1. The persistence R(t) against the rescaled time vari- 
able t/N for different system sizes N = 2™, with n = 12, 15, 
18 and 21 as shown. The thin solid line has a slope of -1.59. 
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FIG. 2. The persistence R{t) for two dimensional lattices 
of sizes N = 2 m x 2 n . The solid straight line has a slope of 
-2.43. The data was averaged over at least 10 separate runs. 
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FIG. 3. R{t) for one dimensional systems of size L = 2 21 
with different interaction rules, to determine the universality 
of the persistence exponent. See the main text for explana- 
tion of the terms "2nd NN", "Half NN" and "Min to 1." The 
solid line has a slope of -1.59. 



FIG. 4. The quantity t{n cov (t))/(N(S(t))) against time for 
1 and 2 dimensional lattices of sizes 2 21 and 2 11 x 2 1() . The 
argument in the main text predicts that this will tend to a 
constant value that is the persistence exponent. For compar- 
ison, horizontal solid lines are plotted for 6 — 1.59 and 2.43, 
which are the exponents for ID and 2D respectively. 




FIG. 5. R(t) for a random nearest neighbour model in 
which the minimum and another K — 1 sites chosen at ran- 
dom from the lattice are assigned new values. The system 
contained N = 2 21 elements, and R(t) was averaged over 100 
runs. The solid lines have slopes of -0.83 and -1.65. 
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FIG. 6. 1 — R(t) versus t/N for supercritical initial condi- 
tions in which every element has a value /, = 1 except for 
a single 'seed' site that has / see d < 1. The ID data was ob- 
tained from a single run on a 2 16 lattice, and the 2D data was 

x 2 10 lattice. The solid lines 



averaged over 10 runs of a 2 
give the known slopes of tall 



0.42 in ID and 0.70 in 2D. 



